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ABSTRACT 



The tracking accuracy of a missile target seeker depends on many variables. For 
a target seeker using a gimbaled platform, an important variable is the friction induced 
bv the preloaded bearings, and by the short wires which connect the target detector 
with the rest of the seeker's electronics. This friction force is nonlinear and sufficiently 
large enough such that accurate position tracking of a target, whether stationary or 
moving, is difficult. Conventional control methods such as P.D. (proportional plus 
derivative) or P.I.D. (proportional plus derivative plus integral) control action can not 
satisfactorily meet the error criteria. To overcome the deficiency of these tw r o methods, 
a model-reference method has been synthesized, relying on idealized predictor corrector 
control to improve the tracking accuracy of the missile seeker. Computer simulations 
using the Dynamic Simulation Language have demonstrated the superior performance 
expected from this method. 
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I. INTRODUCTION 



A. BACKGROUND 

A typical air-to-air missile seeker might consist of a set of focusing optics, a 
detector, signal processing electronics, gyros, rate and position sensors, all mounted on 
a gimbaled platform as illustrated in Figure 1.1. A missile target seeker's function is to 
detect and to track a target, whether stationary or moving, until the target is 
neutralized. There are many variables which affect the accuracy of the seeker. 
Mechanical variables include friction, inertia, location of center of gravity, and 
vibration. Electrical variables include servos, signal processing limitations, and noise. 
Other variables such as thermal noise, external disturbances, electro-magnetic 
interferences, all play an important role in the tracking accuracy of a missile seeker. 
One of the most important considerations is the effect of friction on the tracking 
accuracy. For a seeker with a gimbaled platform, the friction primarily comes from the 
preloaded bearings and from the wires connecting the platform electronics and the 
seeker electronics located off the platform. This type of friction is nonlinear and its 
direction changes in relation to the motion of the target. Excessive friction causes 
"sticking" behavior and impairs accurate tracking. 

Control of the position of a massive object in the presence of friction is not an 
easy task and requires high positional feedback gain with accompanying high drive 
stiffness for accuracy. Not only does this problem arise in target seekers, but also in 
optical tracking telescopes, and robotic manipulators. Classical control techniques for 
accurate tracking have been based on the provision of a drive torque being 
proportional to angular position error with dynamic compensation based on integral 
and, or derivative action. Target position is sensed by the optical apparatus as an 
angular error from the straight line pointing vector. Platform drive torques are then 
commanded to drive the angular pointing error to zero so that the angular position of 
the gimbaled platform automatically aligns with the target. In stabilizing dynamic 
motions, rate feedback from platform is essentially provided by either the gyros or by a 
resolver mounted opposite the torque motor on the platform. The accuracy of tracking 
is limited by the sensitivity of the feedback elements and is further impaired by 
vibration and sensor noise [Ref. 1[ Normally, high positional feedback gain results in a 
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high drive stiffness, but, with the presence of vibration and noise, high stiffness 
increases errors. Modern stabilization techniques are being sought that overcome the 
need for high feedback gains. 

In this thesis, the use of a model based predictor corrector control has been 
demonstrated, including the use of a "stick-slip" friction model for predicting the added 
torque required to compensate for frictional induced inaccuracies. 




B. RECENT DEVELOPMENTS IN CONTROL 

A typical method employed to control the position error of a mechanical system is 
the P.D. (proportional plus derivative) action. This control scheme works well in most 
control applications but it has one drawback, the presence of steady state error. 
Integral action is usually added to eliminate the steady state error, but it is useless in 
the presence of non-linearity elements. Lag lead compensation technique is usually 
employed to overcome non-linearity effects by increasing the gain in the system and 
subsequently lead to increase noise in the system and increase error. 

The use of model-reference controls increased in recent years. It is found that this 
technique lends itself to the problem of controlling complex and often non-linear 
systems. With the advent of modern microprocessor based controls, techniques using 
model following controls [Ref. 2) offer advantages where a large part of the control 
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force and effect can be predicted, rather that relying on feedback methods alone for 
correction of errors. Model Reference Controls have been used as an ideal response 
generator for comparison with actual system response so that plant parameter changes 
can be monitored, identified, and control gains adjusted automatically. Additionally, 
with particular emphasis on nonlinear systems in robotics, the computed torque 
method and sliding mode control design [Ref. 3] have been proposed for providing 
improved model following and tracking system performance. 

C. OBJECTIVE 

The effectiveness of a target tracker depends on its ability to maintain the target 
within its line of sight regardless of the actions taken by the target. This requires a 
high degree of accuracy on the part of the seeker to response to the motion of the 
target. It is demonstrated that friction, among other impedimenta, reduces the ability 
of the seeker to track a target. It is the purpose of this investigation to compare 
tracking performance with both classical P.I.D., P.D., and model-reference control 
actions under the assumption of deterministic signals. The primary disturbance is 
considered to arise from friction which is modeled here to include coulomb friction and 
the elastic effects from wiring harnesses. 

D. METHOD AND APPROACH 

The approach to this problem begins with a study of the friction variables and 
their effect on the tracking error. Then, the performance of the classical control 
actions, P.D., P.I.D. is evaluated. Finally, a model-reference method has been 
synthesized to improve the tracking accuracy of the seeker. All evaluations were done 
by computer simulation using the Dynamic Simulation Language (D.S.L.) developed 
for the I.B.M. 3033 computer. The simulation is limited to a single degree of freedom 
gimbal and it is shown that superior tracking accuracy may be achieved using a model- 
reference system with nonlinear friction force compensation. 
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II. ANALYTICAL MODELING 



A. INTRODUCTION 

This chapter discusses the development of the mathematical model for the target 
seeker dynamics in one degree of freedom and the reference model used for nonlinear 
friction compensation. The analytical modeling begins with the derivation of the 
equations of motion for the platform in the Y-Z plane. The motions in the other 
planes are neglected at this time to reduce the level of complexity. Even with this 
simplification, the information obtained in this study should prove useful in gaining an 
insight to the interaction of friction in a dynamic system. A friction model based on a 
"stick-slip" concept will be developed to model the combined friction effect of the 
bearings and wire harness. This is followed by the presentation of the classical control 
actions. A model-reference method to control the position error will be synthesized to 
improve the tracking accuracy. A brief discussion of the input signal model will 
complete this chapter. 

B. EQUATION OF MOTION 

The platform, together with the optics, detector, gyros, signal processing 
electronics and sensors can be modeled as a simple cylinder with a small degree of out 
of balance. The analysis begins with the determination of the forces acting on the 
system as shown in Figure 2.1. Using Newton's second law of motion and 
D'Alembert's Principle, the equation of motion of the seeker is derived as, 

U m + mr) 9 + F(0.0) + mzr(cos0cos<p — sinBsintp) = T x (eqn 2.1) 

where the torque, T , required to move the platform through an angular displacement, 
0, depends on the position and velocity of the platform and on the motor's 
characteristics. The friction term in this equation is nonlinear. It is a function of 0 
and 0 The terms in z refer to an axial acceleration induced torque arising from some 
small out of balance mass m located at a distance r from the geometric center making 
an angle (p from the Y axis. 
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With a typical d.c. motor, the output torque can be expressed for purposes of 
position control as the sum of proportional error and rate terms, as, 

T x - K,(0 c -e)-K : 0 (eqn 2.2) 

where Kj is the motor gain constant and K, is the velocity gain constant. For an 
input command. 0 c , the required torque can be calculated by Equation 2.2 if the 
displacement 0 and velocity 0 of the platform are known. The damping ratio v ' 2 is 
used to provide proper damping of the system. 




Figure 2.1 Free Body Diagram of the Platform. 



C. FRICTION MODEL 

The main contributors of friction in this platform are the bearings and the 
electrical wires which connect the platform electronics with that located in the support 
structure. A proper friction model must begin with the understanding of how these 
two components behave. The bearing friction is basically rolling friction which is a 
function of the surface roughness of the races and of the angular speed of the platform. 
The wires act as a spring which could be quite stiff when they are bundled, or soft 
when they are lose. Relative motion between wires in a harness causes an additional 
friction force that is dependent on the stiffness of the harness. The combined friction 



from the two components can be modeled as a massless slider with a spring attached to 
it as shown in Figure 2.2. 




Figure 2.2 Friction Model of the Bearings and Wires. 



The force pulling on the spring stretches the spring by a distance 6 until the 
spring s potential force is overcome. Then, the slider begins to move to a distance 0.. 
The friction force then is equal to the product of the spring constant and the net 
motion ( © — © s ) of the spring. The spring constant, k, is a function of the coefficient of 
friction between the slider and the surface, ft, and of the amount of the allowable 
stretch on the spring. 6. A large value for 6 represents a soft spring; and a small value 
for 6 represents a stiff spring. When the direction is reversed, the force pushes on the 
spring, compressing it until the force is greater than the spring's potential force and the 
slider moves in the opposite direction. The friction force of this model is given as. 



F 



M (0-O s ) 
6 



(eqn 2.3) 



When the friction force is less then the friction between the two surfaces, the slider 
does not move until the friction force increases to exceed fl as indicated in 2.4 and 2.5. 

0 s = 0^ F < p feqn 2.4) 
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F 



e. 




- 6 

ABS(F) 



F > n 



(eqn 2.5) 



In reality, a friction force can exists even though there is no relative motion. A closer 
look at equation 2.3 reveals the ability of this model to support a non-zero friction 
force. When 0 $ is zero, that is, the platform is not rotating, the friction force is equal 
to the product of the spring constant and 0. A typical coulomb friction model would 
have a zero friction force when 0 $ is zero. It is the ability of this model to handle 
nonzero friction force which sets it apart from conventional coulomb friction model. 
Once the force exceeds the friction force, the slider starts to move with its direction 
determined by the direction of the force as indicated in Equations 2.5. The response of 
the slider-spring friction model is illustrated in Figure 2.3. 



F 




M - 


' fk > 0 


-6 


/r ■•!/- • " 




1 / * 0 
6 / 


\L 


— •M-' 



Figure 2.3 Response of the Friction Model. 



D. BLOCK DIAGRAM OF BASIC CONTROL SYSTEM 

The behavior of the system can now be put in a more convenient form for control 
analysis and simulation. A block diagram, is constructed from Equations 2.1 and 2.2. 
Figure 2.4 provides an visual interpretation of the control system which allows more 
effective analysis of the problem. 

The command signal, 0 C , is compared with the position feedback signal, 0, to produce 
an error c. This signal is amplified by a feedforward gain K, in the P.D. control 
method. The amplified signal is compared with the velocity feedback signal to produce 
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Figure 2.4 Block Diagram 
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the required torque to drive the platform. However, this torque is reduced by the 
friction torque caused by the bearings and wires. The torque generated by the out of 
balance mass reduces the available torque even further. What is left is a greatly 
reduced torque which is not necessarily high enough to drive the platform dynamics to 
the desired position and a position error exists. P.D. control action simply would not 
produce a zero position error. The shortcoming of this scheme is its lack of history of 
how the error changes with time. The integral action, shown as dotted outlines in 
Figure 2.4, helps to reduce the error by summing errors accumulated up to the most 
recent time, then an amplified average of this error is added to produce the required 
torque. The overall effect is better error reporting and correcting. Ideally, this method 
would produce zero position error. Both the P.D. and P.I.D. control suffer noise 
sensitivity problems. When a noise, either internal or external, is introduced to the 
system it is amplified by the control actions and accurate tracking becomes difficult. 

E. MODEL-REFERENCE CONTROL 

Consider an ideal system which has the same basic inertial characteristics as the 
platform. This is refered to as a model. Assuming the model is perfectly balanced and 
neglecting the effect of friction, the ideal model dynamics is expressed by Equation 2.6. 



If this model is rotated by a motor with the output torque expressed as the sum of a 
positional error and rate feedback, the torque equation is given as, 



where U is the ideal motor output torque, 0 c is the command signal as defined before. 
X and X denotes the position and velocity of the model, respectively. Gj is the 
positional error gain and G 2 is the velocity feedback gain. The value of G ( is limited 
only by the deliverable torque from the torque motor and G, follows the same 
damping rule as K 2 but without the problem of signal noise generation. The model 
friction due to the bearings and wires can be modeled as before. Using the notation X 
in place of 0, X in place of 0, n m and 6 m replace n and 6. the model friction is 
expressed by Equation 2.8. 



U 



X 




(eqn 2.6) 



U = G,(0 c -X)-G 2 X 



(eqn 2.7) 
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(eqn 2.8) 



u (X-x ) 

p — "m s 

m § 

m 



X s = X s 



X = X - 8 22 

5 s m ABS{F m ) 



F m < ^ n 2 ' 9 ) 

F m > M m (eqn 2.10) 



In order to overcome the friction torque and drive the model to 0 c , the torque needed 
must be the sum of the two torques described above. The predicted torque is given 
below. 



T xpred = L ' + F m 



(eqn 2.1 1) 



To correct for any positional and velocity errors that might exist between the model 
and the platform, an additional correction torque based on the amplified error signals 
are added to the predicted torque to yield the total required torque, expressed in 
Equation 2.12 as, 

T x = T xpred + K l( X " 0 ) + K 2 (X-0) (eqn 2.12) 



where the second term corrects the positional error and the third term corrects the 
velocity error. 

This software based model is described in Figure 2.5. It is assumed that the torque 
to drive the platform in following command signals can be expressed as the sum of the 
torque to drive an ideal platform with no friction plus the torque required to overcome 
friction as established by the ideal platform model. Since the actual platform motion 
may differ slightly from that of the ideal, a correction torque is added based on a 
proportional and derivative action applied to motion errors. 

F. INPUT SIGNAL MODEL 

The response of the system to a step input and a sine input is simulated using the 
D.S.L. simulation program. An input command is 0.2 radians for all simulations. This 
is equivalent to commanding the seeker to move 11.459 degrees. The step input is 
roughly equal to an initial target acquisition phase which the seeker is required to lock 
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D 




Figure 2.5 Model- Reference Control Action 
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on to the target or when the target is stationary. The sinusoidal input simulates the 
maneuvering of the target to avoid lock on. The target is assumed to be moving at 2.5 
Hz. This type of maneuver is extreme because most targets are not able to evade this 
quickly. However, most missiles flex and rotate while moving through the air, which 
makes this assumption more realistic. 



III. COMPUTER SIMULATION METHOD 



A. INTRODUCTION 

Dynamic Simulation Language (D.S.L.) is a FORTRAN-based simulation 
language for simulation of continuous systems. It's strength lies in its built in 
functions which allow the composition and simulation of any physical systems. Some 
of the built in functions include integrators, function generators, non-linear functions, 
probability distribution functions, and linear transfer functions; allow easy construction 
and simulation of the system without heavy programing. A more comprehensive look 
at the capabilities of D.S.L. can be found in References 4 and 5. 

B. GEOMETRIC DESCRIPTION 

The platform was modeled as an cylinder two inches in diameter and three and a 
half inches in height. The total mass of the cylinder was 0.0029 slugs plus a 10% out 
of balance mass (0.0003 slugs). The inertia of the cylinder was 0.0037 lb-in-s 2 The 
initial distance of the out of balance mass from the geometric center was taken to be 
0.1 inches in the positive V direction and 0.2 inches in the positive Z direction. This 
distance as well as the friction and control parameters were varied to determine their 
effects on the performance of the seeker. 

C. ACCELERATION SIMULATION 

The effect of acceleration is simulated under an assumed flight profile. The 
maximum acceleration subjected to the platform is assumed to be an exponential 
function describes by Equation 3.1 and 3.2. 



Z = Z m (l-EXP{- t/tj) 


t < r, 


(eqn 3.1) 


Z = Z m EXP{- t;t 2 ) 


t> 4t, 


(eqn 3.2) 



where and t 2 are some arbitrary constants. 

D. DSL CODING 

The flow chart shown in Figure 3.1 illustrates the basic structure of the D.S.L. 
program used in the simulation of the seeker. The program is divided into segments: 
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TITLE, CONSTANT, PARAMETER, INITIAL, DYNAMIC, DERIVATIVE. AND 
TERMINAL. The TITLE segment named the program. The CONST segment defines 
the constants in the program which include the mass properties of the seeker and the 
maximum acceleration. PARAM segment defines the constants which will be changed 
from one simulation run to the next. It defines the friction characteristics of the seeker 
model, input signals, time constants, gain constants, and the geometry of the out of 
balance mass. The INITIAL segment initializes the variables and calculates the values 
of the out of balance geometry and velocity feedback constant. The DYNAMIC 
segment calculates the acceleration profile of the seeker and the input signals. This 
segment is computed at each time steps. The DERIVATIVE segment is the main body 
of the program. It consists of the description of the seeker dynamics as well as the 
control actions being used. Finally, the last segment, TERMINAL, contains the 
commands with regard to the total simulation time desired as well as the printing and 
plotting information. 

Appendix A presents the DSL programming codes for the simulation of the 
classical P.D. and P.I.D. control actions and Appendix B presents the DSL 
programming codes for the model-reference control method. The symbol "ff" means 
the line is not used for current simulation. 
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Figure 3.1 D.S.L. Program Flowchart 



IV. SIMULATION RESULTS 



A. INTRODUCTION 

This chapter presents the results of the simulation. The effect of friction, control 
parameters, and control schemes on the tracking accuracy are illustrated by 
accompanying figures. Except as indicated, all figures shown the command signal, 0 c , 
as a series of short dashes and the actual response, 0 (TH). as a solid line. The actual 
torque, T (TX) was indicated as medium length dash lines. Most of the figures reflect 
the response to the sinusoidal input, and also contain the system error (THE) as 
indicated by long dash lines. Other information in the figures include the error 
between the response and the slider (THERR). This piece of information indicates 
how the slider from the friction model relates to the response. The pound(lb), inch(in), 
second(sec) system was used through out the simulation process. The unit of torque 
was lb-in. The unit of angular displacement was in radians (rad). The velocity and 
acceleration were in rad sec and rad sec 2 , respectively. Mass was expressed in lb- 
sec 2 in. The command signal amplitude was 0.2 radians (11.459 °) and the simulation 
time lasted one second. 

B. EFFECT OF FRICTION ON ACCURACY 

This section discusses the effect of friction on the system. The system with no 
friction was simulated and its result was compared with the same system with various 
friction level. In order to provide a meaningful result, all the simulation variables 
except ft must remain constant. The control variables K L and 5 were set at 15 lb- 
in rad and 0.05 in, respectively. The value of K, was set at 0.5 lb-in-sec, 'rad. Figures 
4.1 through 4.8 are related to this subject. The response of an ideal, balanced, system 
with no friction, to a step input is shown in Figure 4.1 The rise time 1 of the ideal 
response was 0.12 seconds with no steady state error. Using the P.D. control scheme 
with 4 equals to 0.1, 0.2, and 0.3; the steady state error increased from 0.00667 to 
0.01333 and to 0.02 radians, respectively. The performance of the P.l.D. control action 
was better. The worst steady state error occured when n equal to 0.3 lb-in was only 



l The rise time w r as taken as the time required for the system to achieve 95% of 

the steady state value. 
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0.00424 radians. However, this better accuracy was at the expense of rise time which 
was 0.24 second. This was twice as slow as in the ideal case. A comparison of 
position error for different friction level at different gain Kj was shown in Figure 4.8. 
This figure also shown some interesting information. It was expected the position 
error would increase with increase friction at any given set of circumstances. However, 
the figure indicated that although it was true in general there were cases when an 
increased did not produce a higher error. For example, when equal to 10 lb- 
in rad the position error for H equal to 0.2 lb-in had a higher position error (2.5 milli- 
radians) compared to 0.84 milli-radians for n equal to 0.3 lb-in. This was due to the 
nonlinearity of the system. A particular K 2 combined with a suitable Kj and 5 would 
produce a better error. 
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Figure 4.1 Balanced System, PD Control. No Friction 
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Figure 4.2 H = 0.1. Kj = 15 0. K, = 0.5. 6 = 0.U5 
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Figure 4.3 n = 0.2, K,= 15.0, K 2 = 0.5, 6 = 0.05 
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figure 4.4 p = 0.3, Kj - 15.0. K, = 05. 6 = 0.05 
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Figure 4.5 fi = 0.1. K, = 15.0, K 2 «0.5. 5 = 0.05, PID Control Tj« 1.0 
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Figure 4.6 n = 0.2. K, = 15.0, K 2 = 0.5, 5 = 0.05, PID Control T= 1.0 
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Figure 4.7 n = 0.3. K, = 15.0. K, = 0.5. 6 = 0.05, P1D Control T= 1.0 
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Figure 4.S Position Error Vs Gain Constant K ( 
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C. EFFECT OF THE K , 

The usual method of minimizing the influence of disturbance forces such as the 
frictional load is to increase the feedback gain, K 1? so that system "stiffness" is 
increased. However, while this works with clean signals, the effects of sensor noise 
provides a limit on the useful range of Kj. The deliverable torque by the motor also 
limits the values of Kj. Figures 4.9 through 4.14 provide graphical results to different 
K^s. A high value drastically reduces the position error of the system as shown in 
Figure 4.8. The steady state error decreased from 73.88% to 10% when Kj increased 
from 2 lb-in rad to 15 lb-in rad. Like wise, the positional error to a sinusoidal input 
decreased from 0.12559 radians to 0.01764 radians when Kj increased from 15 lb-in rad 
to 1000 lb-in rad. The reponse followed the command signal rather well at the higher 
Kj value with a phase shift of less than 0.01 seconds. The torque delivered to the 
system was also higher, an increase of 65%, compared with the lower Kj value of 15 
lb-in rad. 

D. EFFECT OF SPRING CONSTANT 

Figures 4.15 through 4.20 shown the effect of spring stiffness to position error. 
The spring constant, 5, appeared to have little effect on the steady state error but it did 
affect the torque delivery 7 to the platform. A very stiff spring, small 5, required more 
torque than a softer spring as expected. Indeed, it was shown in Figures 4.15 and 4.16 
However, this logic did not work in Figure 4.17. The torque delivery in this case was 
higher. The effect on a sinusoidal input was a little different. The stiffer spring 
performed poorer than the softer spring. The maximum error was 0.18527 radians for 
the stiff spring and 0.14129 radians for the softer spring. However, the softer spring 
followed the signal a little better than the stiff spring. The time lag was 0.07 second for 
the stiff spring and 0.04 second for the soft spring. Compares with the ideal case where 
there was no friction, the maximum error was 0.12092 radians with a time lag of 0.04 
seconds as shown in Figure 4.18. 
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Figure 4.9 n = 0. 3. K[ = 2 . 0 . 5 = 0.05 
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figure 4.10 n = 0.3, K, = 5.0, 5 = 0.05 
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Figure 4.11 n = 0.3, K, = 15.0. d = 0.05 
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figure 4.12 H = 0.3, Kj = 15.0. 0 = 0.05, Sin Input 
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Figure 4.14 H = 0.3, K, = 1000.0, 6 = 0.05, Sin Input 
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Figure 4.15 (1 = 0.3. Kj = 
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Figure*!. 16 fl = 0.3. K ( = 5.0,6 = 0.01 
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E. EFFECT OF INTEGRAL GAIN TIME CONSTANT 

The benefit of an integral action was demonstrated in Section B. This section 
studied the effect of the integral time constant. When the time constant was set to 
unity (Figure 4.5), the response overshot by 0.00188 radians (0.94%) and then settled 
back to 0.20104 radians at the end of one second. With the time constant increases to 
2 and 3 seconds, the response did not exhibit any overshoot as shown in Figures 4.21 
and 4.22. It did indicate that as the time constant increases, the steady state error also 
increases, even though the increase was small. 

F. EFFECT OF G , 

This section as well as the next section focuses on the control variables Gj and 
prediction on the accuracy of the seeker using the model-reference technique. It was 
shown in previous sections that the sinusoidal input produced an error due to the 
inability of the seeker to track accurately, even in the case where friction was not 
present. This was mainly due to the inability of the P.D. or P.I.D. control action's 
failure to correct the inertia effect of the platform in a timely manner. The effect of Gj 
was illustrated in Figures 4.23 through 4.25. In Figure 4.23, the maximum error was 
0.09086 radians and a time lag of 0.03 seconds with Kj and G t both equal to 30. It is 
interesting to observe that the performance did not degrade even though the value of 
Kj was reduced by half provided the value of Gj was high. In fact, the performance 
can be improved by increasing G, As shown in Figure 4.24 and 4.25, both the 
maximum error and time lag were drastically reduced to 0.02417 rad and 0.01 rad as Gj 
increases to 500 and 3000, respectively. The amount of time lag also cut by a third to 
less than 0.01 second. Although the value of Gj was arbitrary, it was not without 
bound. The maximum usable value depends on the size of the torque motor. 
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Figure 4.17 n = 0.3 f K.- 5.0, 0 = 0.1 
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Balanced System. PD Control, No Friction. Sin Input 
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Figure 4. IS 
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Figure 4.19 fi = 0.3, K 



5.0, 6-0.01, Sin Input 
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M = 0.3, K, = 5.0. 6 = 0.1. Sin Input 
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Figure 4.21 jl = 0.1. Kj = 15.0, 5 = 0.05. FID Control, 12 = 2 
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Figure 4.22 H = 0.1, Kj = 15.0, 6 = 0.05, PID Control, T= 3 
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Figure 4.23 n = 0.3. K 



30.0. <> = 0.05, ft = 0.3. 5 = 0.05. G. = 30.0 
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Figure *1.24 n = 0.3. K t = 30.0, 6 = 0.05, M m = 0.3, 6 m = 0.05, G, = 500.0 



G. EFFECT OF \i PREDICTION ERROR 

The main advantage of the model-reference control scheme is its ability to predict 
the response of the system using an ideal model. Figures 4.26 through 4.29 study the 
sensitivity to an erroneous prediction of the friction level in the system. Friction levels 
of p equal to 0.1 to 0.5 with an increment of 0.1 were simulated. The results indicated 
that a wrong prediction did make a difference. The maximum errors in order of 
increasing p were 0.01692, 0.01862, 0.01, 0.01368, and 0.01845. A smaller system 
friction than predicted tends to cause overshoot while a larger system friction causes 
the response to be a little short. In the sinusoidal input response, the corrective torque 
was trying to correct the overshoot but it's effectiveness was quite limited. This 
problem can be corrected by adaptive control action. Very often, the model-reference 
action was used with adaptive control. 

H, EFFECT OF UNBALANCED SYSTEM 

The effect of the out-of-balance appeared to be quite minimal as shown in Figures 
4.30 through 4.32. The response typically reached a maximum values early on, usually 
within the first 0.12 seconds, then it droped slowly as the acceleration reached 
maximum acceleration. After the platform reached the maximum acceleration and 
started to decelerate, the response rose again. As the distance between the center of 
gravity and the out of balance mass increased, the accuracy decreased. Therefore. It is 
important to consider the effect of out-of-balance mass in platform design. 
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Figure 4.25 
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Figure 4.27 *1 = 0.2. K,= 15.0, 6 = 0.05, M m = 0.3, 6 m = 0.05, Gj = 3000.0 
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Figure 4.2S u = 0.4. K, = 15.0. 6 = 0.05, H =0.3.6 = o.u5. G, = 3000.0 
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figure 4.29 



= 0.5, = 15.0, 6 = 0.u5. }i |n = 0.3, S m = 0.05. G, = 3000.0 
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Figure 4.30 fi = 0.3, K, = 15.0, 5 = 0.05, PD Control. Z m = 3S64.0, Y o = 0.1, Z o = 0.2 
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Figure 4.31 fi = 0.3, K,= 15.0, 5 = 0.05, PD Control. Z m = 3864.0, y o = 0.2. Z o = 0.4 
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Figure 4.32 M = 0.3, Kj= 15.0, 6 = 0.05, PD Control, Z m = 3864.0, y o = 0.2, Z o = 0.6 
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I. EFFECT OF ACCELERATION 

The effect of acceleration was similar to that of the C.G. location. The response 
began to drop as the acceleration approached maximum and recovered when the 
platform started to decelerate. This phenomenon was true for both the P.D. and 
model-reference control schemes shown in Figures 4.33 through 4.36 with either step 
input or sinusoidal input. However, the model-reference action allowed a much 
quicker recovery than the P.D. action. This was due to the correction torque provided 
by the model. Furthermore, the P.D. action had a much greater positional error while 
responsed to a step input and a slower response to a sinusoidal input. 
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Figure 4.33 p = 0.3, K t = 15.0, 6 = 0.05. Z m = 3864.0, PD Control, step response 
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Figure 4.34 



p = 0.3, Kj = 15.0, 6 = 0.05, Z (n = 3S64.0, FD Control, sin response 
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Figure 4.35 n = 0.3. K, = 15.0. 6 = 0.05, Z m = 3864.0. M-F Control, step response 
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Figure 4.36 n = 0.3, K t = 15.0, 6 = 0.05. Z m = 3864.0, M-F Control, sin response 
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V. DISCUSSION AND RECOMMENDATION 



A. INTRODUCTION 

This chapter summaries the results obtained from the D.S.L. simulation and make 
recommendations for future work in this area. 

B. DISCUSSION 

The model-reference control method described here offered better tracking 
accuracy than either the P.D. or P.I.D. actions. It also minimized the problem of 
noisy feedback signals. This software based model utilizes an idealized model as a 
reference to predict the torque needed to overcome friction and gimbal dynamics and 
produce zero position error. First, the amount of torque needed to drive an idealized 
platform with no friction is predicted. Then, the torque required to overcome friction 
is determined from the friction model described in Chapter II. The sum of these two 
produced a predicted torque T ed . In addition, the idealized platform velocity X is 
compared with the actual velocity 0 and the error is amplified. Similarly, the ideal 
position X is compared with the actual position 0 and the error is amplified also. A 
corrective action occurs when these amplified error signals are added to T xpred to yield 
the total required torque T . The result is a highly tuned, high accuracy, control 
action. Since the reference model is ideal, it does not have the draw back of large noise 
amplification as the other two control schemes do. Therefore, it is reasonable to 
expect better target tracking, even if the noise level is high. 

The simulation results clearly shown the superior tracking accuracy provided by 
the model-reference method. The level of accuracy achieved is adequate for most 
target tracking seekers. Although this study was limited to the Y-Z plane only, this 
exercise did provide a better and clearer understanding of how friction and control 
variables affect the performance of the seeker. Of course, the performance of a missile 
seeker is affected by many other variables, not discussed here. The effect of vibration 
as the missile moves through the air, or the effect of air blast in the vicinity of the 
missile are examples of disturbance which will play an important role in the tracking 
accuracy of the missile. The friction model presented here is more realistic than simple 
coulomb friction model which produces an underfined friction force when 0 $ is equal to 
zero. The model presented here does allow non-zero friction force under stationary' 
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conditions. Improvements can be made on the model-reference method to have 
adaptive capability to take into consideration the external disturbances, and 
variabilities in actual friction. 

C. RECOMMENDATION 

It is recommended that an experimental follow up be made to compare with the 
results obtained by analytical means. Frictional effect is not the only element affect 
seeker performance. There are other problem areas which also deserve the attention 
given here; among them vibration, shock, and random noise effects. It will be 
interesting to subject the model-reference control method described here to these 
external disturbances to see how well or how poor the current method handles these 
problem. These disturbances can be incorporated into the model-reference algorithm 
so their effects will be minimal. There are also an abundance of research in adaptive 
control. This type of control action enables a control system such as the one described 
here the ability to "learn and adjust" the model to different signal input. This latest 
technique is most versatile and would definitely contribute to the tracking accuracy of 
the seeker. 
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APPENDIX A 

DSL LISTINGS (PD, PID CONTROL) 



TITLE GIMBAL DYNAMICS 
* 

* * 

* AUTHOR: JOSEPH CHAN * 

* ABSTRACT: THIS PROGRAM SIMULATES THE MOTION OF A SINGLE * 

* GIMBAL SYSTEM WHILE SUBJECTS TO STEP AND SINUSOIDAL * 

* EXCITATIONS. PD AND PID CONTROL ACTION. * 

* * 

* 



* 

* 

SYMBOLS : * 

* 

BM : MASS OF THE CYLINDER. * 

DELTA : CHAR. DISP FOR COULOMB FRICTION MODEL. * 

F : BEARING FRICTION FORCE. * 

IM : MASS MOMENT OF INERTIA OF CYLINDER. * 

K1 : TORQUE GAIN CONSTANT. * 

K2 : VELOCITY FEEDBACK CONSTANT. * 

MU : COEFFICIENT OF FRICTION. * 

FERIOD : PERIOD OF THE SINE INPUT. * 

PH : ANGLE BETWEEN THE Y-AXIS AND R. * 

R : DISTANCE TO THE OFF-CENTER MASS FROM THE ORIGIN. * 

SM : OUT OF BALANCE MASS. * 

TAU1 : ACCELERATION TIME CONSTANT (INITIAL PHASE). * 

TAU2 : ACCELERATION TIME CONSTANT (FINAL PHASE). * 

TH : ANGULAR DISPLACEMENT OF THE CYLINDER. * 

THO : INITIAL ANGLE ROTATION. * 

TH2DOT : ANGULAR ACCELERATION. * 

THCOM : COMMAND ANGLE INPUT. * 

THDO : INITIAL ANGULAR VELOCITY. * 

THDOT : ANGULAR VELOCITY. * 

THE : ANGULAR DISPLACEMENT OF THE CYLINDER. * 

THS : POSITION OF THE SLIDER. (FRICTION MODEL) * 

TIME : SIMULATION TIME. * 

TX : MOTOR TORQUE. * 

YO : DISTANCE TO OFF CENTER-MASS IN THE Y DIRECTION. * 

ZO : DISTANCE TO OFF CENTER-MASS IN THE Z DIRECTION. * 

Z2DOT : LINEAR ACCELERATION. * 

ZM : MAXIMUM LINEAR ACCELERATION. * 

X * 

********************************************************************** 

* GIMBAL PARAMETERS 
********************************************************************* 

* 

CONST BM=0. 00290, SM=0. 00030, ZM=0.0000, IM=0. 00370 

* LB-S**2/IN LB-S**2/IN IN/S**2 LB-IN-S**2 



* 

* 

* 

X 

* 

* 

* 

* 

* 

X 

* 

X 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

X 

* 

* 

X 

* 

* 

* 



PARAM Kl=15.000, 
LB-IN/RAD 



MU=0 . 3 , 



DELTA=0 . 05 , 
IN 



PARAM Y0=0 . 

* IN 

* 

INITIAL 

R 

PH 

K2 

F 

TH 

THO 



20=0.2, TAU1=0 . 5 , TAU2=0 . 5 , 
IN SEC SEC 



= SQRT (Y0**2+Z0**2 ) 
= ACOS (YO/R) 

= SQRT(8*IM*K1) 

= 0.0 
= 0.0 
= 0.0 



PERIOD=0 . 4 
SEC/CYC 

THMAX=0 . 2 
RAD 
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* * 



THDO 


= 


0.0 


THDOT 


= 


0.0 


THE 


= 


0.0 


THEO 


= 


0.0 


THERR 


= 


0.0 


THS 


= 


0.0 



DYNAMIC 

* 

•k k ***** ************** * ** ****************** **** ****** ** ******* ********* 

* ACCELERATION PROFILE * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k 

T1 = TIME-4. 0*TAU1 

Z2DOT = ZM*(1. 0-EXP (-TIME/TAU1 ) ) 

IF (TIME .GT . (4 .0*TAU1 ) ) THEN 
Z2DOT = ZM*EXP(-T1/TAU2) 

END IF 



* 

* 

* 

# 

k 

k 

k 



SINUSOIDAL EXCITATION 

THCOM=THMAX*SIN(2*PI*TIME/PERIOD) 
STEP INPUT 

THCOM=THMAX 



DERIVATIVE 

NOSORT 

k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* BEARING FRICTION MODEL * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k 

THERR = TH-THS 

F = MU* (THERR) /DELTA 

IF ( (ABS(F) ) . LT.MU) THS = THS 

IF ( (ABS(F) ) .GT.MU) THS = TH-DELTA*F / ABS ( F ) 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* MOTOR TORQUE MODEL * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k 

THE = THCOM-TH 



* 

* 

* 

# 

* 

* 

* 



PD FEEDBACK 
TX 



= K1*(THE) -K2*(THDOT) 



PID FEEDBACK 



TX = Kl*(THE)-K2*THDOT+(Kl/TI )*Y 

TH2DOT = (TX-F-SM*Z2DOT*R*(COS (TH)*COS(PH) -SIN(TH) ... 

*SIN(PH) ) )/(lM+SM*R**2) 

THDOT = INTGRLiTHDO ,TH2DOT) 

TH = INTGRL(THO, THDOT) 

* Y = INTGRL (THEO , THE ) 

METHOD RK5 

CONTRL FINTIM=1 .00 , DELT=0.001 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* PRINT RESULTS * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k 

SAVE 0.001, TX, F, THCOM , TH, THE, THERR 
PRINT 0.010, TX, F, THCOM, TH, THE, THERR 

k 

* THETA VS TIME 

* 

#GRAPH (A,DE=TEK618) TIME (LE=8 , UNIT= ' SEC 1 ) TH(LO=- . 3 , SC= . 1 ,NI=6 , ... 
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# UN= 1 RAD ' ) THCOM( AX=OMIT , L0=- .3,SC=0.1 , NI=6 , LI=4 . UN= 1 RAD 1 ) .. 

# THE ( L0=- . 3 , SC= . 1 ,NI=6,LI=2) , TX(LI=3 , UN= 1 LB- IN ' ) 

GRAPH (A,DE=TEK618) TIME(LE=8 ,UNIT= ' SEC ' ) TH(LO=0 . 0 , SC= . 04 ,NI=3 , . 
UN= ' RAD ' ) THCOM(AX=OMIT,LO=0.0,SC=0.04,NI=8,LI=4,UN='RAD' ) 
TX(LI=3 , UN= 1 LB" IN ' ) 

LABEL (A) MU=0 . 3 , Kl=15.0, DELTA=0.05 

# THETA VS TIME WITH THETA COMMAND AS REFERENCE (SINE INPUT) 

#GRAPH (A,DE=TEK618) TIME (LE=8 , UNIT= ' SEC 1 ) TH(LO=-0 . 3 , SC=0 . 1 ,NI=6 , 

# UN= ' RAD 1 ) THCOM(LO=-0 .3 ,SC=0 . 1 ,NI=6 , LI=4 ,UN= 1 RAD 1 ) 

X 

# HYSTERISIS 

X 

#GRAPH (B,DE=TEK618) TH(LE=8,P0=2,4,AX=LIN,DRAW,UN=RAD) F(UN='LB') 
#LABEL (A,B) MU=0 . 25 , Kl=2.00, K2=0.05, DELTA=0.050 

# TIME RESPONSE 

•k 

#GRAPH (C , DE=TEK618) TIME (LE=8 ,UN= 1 SEC ' ) TH(LO=- . 3 , SC= . 1 .NI=6 , ... 

# LI = 1 ,UN= ' RAD 1 ) THS (LO=- .3 , SC= . 1 ,NI=6 , LI=2 ,UN= ' RAD ' ) ... 

# F (LO=- . 3 , SC= . 1 ,NI=6 , LI=4 , UN= 1 LB 1 ) 

#LABEL (C) TIME HISTORY 

X 

END 

STOP 
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APPENDIX B 

DSL LISTINGS (PC CONTROL) 



TITLE GIMBAL DYNAMICS2 

k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk^ 
* * 

* AUTHOR: JOSEPH CHAN * 

* ABSTRACT: THIS PROGRAM SIMULATES THE MOTION OF A SINGLE * 

* GIMBAL SYSTEM SUBJECTS TO STEP AND SINUSOIDAL EXCITATIONS * 

* USING PREDICTION CORRECTION CONTROL ACTION. * 

k k 



7^ — — . 

k 




k 


SYMBOLS : 


k 




k 


BM 


k 


DELTA 


k 


DELTAM 


k 


F 


k 


FM 


k 


G1 


k 


G2 


k 


IM 


k 


K1 


k 


K2 


k 


MU 


k 


MUM 


k 


PH 


k 


R 


k 


SM 


k 


TAU1 


k 


TAU2 


k 


TH 


k 


THO 


k 


TH2DOT 


k 


THCOM 


k 


THDO 


k 


THDOT 


k 


THE 


k 


THS 


k 


TIME 


k 


TX 


k 


TXPRED 


k 


X 


k 


XO 


k 


XI 


k 


X2D0T 


k 


XDOT 


k 


XDOTO 


k 


xs 


k 


YO 


k 


ZO 


k 


Z2DOT 


k 


ZM 



MASS OF THE CYLINDER. 

CHAR. DISP FOR COULOMB FRICTION MODEL. 

CHAR. DISP FOR COULOMB FRICTION (REFERENCE MODEL). 
BEARING FRICTION FORCE. 

BEARING FRICTION FORCE (REFERENCE MODEL). 

REFERENCE MODEL TRANSFER FUNCTION. 

REFERENCE MODEL TRANSFER FUNCTION. 

MASS MOMENT OF INERTIA OF CYLINDER. 

TORQUE GAIN CONSTANT. 

VELOCITY FEEDBACK CONSTANT. 

COEFFICIENT OF FRICTION. 

COEFFICIENT OF FRICTION (REFERENCE MODEL). 

ANGLE BETWEEN THE Y-AXIS AND R. 

DISTANCE TO THE OFF-CENTER MASS FROM THE ORIGIN. 
OUT OF BALANCE MASS. 

ACCELERATION TIME CONSTANT (INITIAL PHASE). 
ACCELERATION TIME CONSTANT (FINAL PHASE). 

ANGULAR DISPLACEMENT OF THE CYLINDER. 

INITIAL ANGLE ROTATION. 

ANGULAR ACCELERATION. 

CONTROL ANGLE INPUT. 

INITIAL ANGULAR VELOCITY. 

ANGULAR VELOCITY. 

ANGULAR DISPLACEMENT OF THE CYLINDER. 

POSITION OF THE SLIDER. (FRICTION MODEL) 
SIMULATION TIME. 

MOTOR TORQUE. 

MOTOR TORQUE (PREDICTED) 

DISTANCE OF MASS AT ANY TIME T. (FRICTION MODEL) 
INITIAL DISTANCE OF REFERENCE MODEL. 

DISTANCE MOVED BY THE SLIDER. (FRICTION MODEL) 
ACCELERATION OF REFERENCE MODEL. 

VELOCITY OF REFERENCE MODEL. 

INITIAL VELOCITY OF REFERENCE MODEL. 

POSITION OF SLIDER (REFERENCE MODEL). 

DISTANCE TO OFF CENTER-MASS IN THE Y DIRECTION. 
DISTANCE TO OFF CENTER-MASS IN THE Z DIRECTION. 
LINEAR ACCELERATION. 

MAXIMUM LINEAR ACCELERATION. 



k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 

k 



kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

k GIMBAL PARAMETERS * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkr:kk 



k 



CONST BM=0. 00290, 
* LB-S**2/ IN 



SM=0. 00030, ZM=0 . 0000 , IM=0. 00370 

LB-S**2/ IN IN/S**2 LB-IN-S**2 



PARAM K1=1500Q0 , MU=0.3, DELTA=0.05 



PERI0D=0 .4 
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LB- IN/ RAD 
PARAM Gl=500 . 00 , 

X 

k 

PARAM Y0=0.1, 

* IN 



MUM=0 . 3 , 



INITIAL 



IN 

DELTAM=0 .05 , 
IN 



SEC/CYC 

THMAX=0 . 2 
RAD 



20=0.2, TAU1=0 . 5 , TAU2=0 . 5 
IN SEC SEC 



R = SQRT(Y0**2+Z0**2) 

PH = ACOS ( YO/R) 

K2 = SQRT(8*IM*K1) 

G2 = SQRT(8*IM*G1) 

F =0.0 

TH =0.0 
THO =0.0 
THDO =0.0 
THDOT =0.0 
THE =0.0 
THEO =0.0 
THERR =0.0 
THS =0.0 
X =0.0 

XO =0.0 
XDOT =0.0 
XDOTO =0.0 
X2D0T =0.0 
XS =0.0 



DYNAMIC 

X 

* x * * * * x * * * * * * * * * * * * k k k * * * * * * A A * * * * k * * * * * * * x * x * * x * & k * * * * * A * * k A A * * * x * x * k 

* ACCELERATION PROFILE * 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



T1 = TIME-4. 0*TAU1 

Z2DCT = ZM* (1 . 0-EXP ( -TIME/TAU1 ) ) 

IF (TIKE.GT. (4.0*TAU1)) THEN 
Z2DOT = ZM*EXP( -T1/TAU2 ) 

ENDIF 

k 

* SINUSOIDAL EXCITATION 

X 

* THCOM=THMAX*SIN(2*PI*TIME/PERIOD) 

X 

* STEP INPUT 

k 

THCOM=THMAX 

k 

DERIVATIVE 

NOSCRT 

X 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* MODEL REFERENCE DYNAMICS 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



U = (THCOM-X) *G1-G2*XD0T 

X2DOT = U/(IM+SM*R**2) 

FM = MUM/DELTAM* ( X-XS ) 

IF ( (ABS(FM)) .LT .MUM) XS = XS 

IF ( (ABS (FM) ) . GT .MUM) XS = X-DELTAM*FM/ABS (FM) 

TXPRED = U+FM 

k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* BEARING FRICTION MODEL 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



THERR = TH-THS 

F = MU* (THERR) /DELTA 

IF ( (ABS (F) ) .LT .MU) THS = THS 
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X- X- X- X- X- 



^ IF ( (ABS(F) ) .GT .MU) THS = TH-DELTA*F/ABS(F) 

********************************************************************** 

* MOTOR TORQUE MODEL 

* (PREDICTOR CORRECTOR METHOD) 

************* ************** ******************************************* 

* 

THE = THCOM-TH 

TX = TXPRED+K1* ( X-TH ) +K2* ( XDOT-THDOT ) 

IF (TX.GT.55.0) TX=55.0 

TH2DOT = (TX-F-SM*Z2DOT*R*(COS(TH)*COS(PH)-SIN(TH) ... 

*SIN(PH) ) )/ (IM+SM*R**2) 

THDOT = INTGRL ( THDO , TH2DOT ) 

TH = INTGRL THO, THDOT) 

X = INTGRL (XO,XDOT) 

XDOT = INTGRL (XDOTO ,X2DOT) 

METHOD RK5 

CONTRL FINTIM=1 . 00 , DELT=0.0001 

********************************************************************** 

* PRINT RESULTS * 

********************************************************************** 

* 

SAVE 0.001, TX, TXPRED , F, THCOH, TH , THE, THERR 
PRINT 0.010, TX, TXPRED, F, THCOM, TH, THE, THERR 

* THETA VS TIME 

* 

#GRAPH (A,DE=TEK618) TIME (LE=8 , UNIT= 1 SEC 1 ) TH(LO=- . 3 , SC= . 1 ,NI=6 , ... 

# UN= 1 RAD ' ) THCOM (AX=OMIT , LO=- . 3 , SC=0 . 1 , NI=6 ,LI=4 , UN= 1 RAD 1 ) ... 

# THERR (LO=- . 3 . SC=. 1 ,NI=6 , LI=2 ) , TX(LI=3) 

GRAPH (A,DE=TEK618) TIME (LE=8 , UNIT= ' SEC ' ) TH(LO=0 .0 , SC= .04 ,NI=8 , ... 

UN= 1 RAD 1 ) THCOM (AX=OMIT , L0=0 . 0 , SC=0 . 04 ,NI=8 , LI=4 , UN= 1 RAD ' ) ... 
TX(LI=3) 

LABEL (A) MU=0 . 3 , Kl=1000, DELTA=0.05, MUM=0.3, DELTAM=0 . 05 , Gl=1000 

* THETA VS TIME WITH THETA COMMAND AS REFERENCE (SINE INPUT) 

#GRAPH (A,DE=TEK618) TIME (LE=8 , UNIT= ' SEC 1 ) TH(LO=-0 . 3 , SC=0 . 1 ,NI=6 , ... 

# UN= 1 RAD 1 ) THCOM (L0=-0 . 3 , SC=0 . 1 ,NI=6 ,LI=4 , UN= 1 RAD 1 ) 

* HYSTERISIS 

* 

ttGRAPH (B ,DE=TEK618) TH(LE=8,PO=2,4,AX=LIN,DRAW,UN=RAD) F(UN='LB') 

# LABEL (A,B) MU=0.25, Kl=2.00, K2=0.05, DELTA=0.050 

* TIME HISTORY 

* 

#GRAPH (C , DE=TEK618 ) TIME (LE=8 , UN= ' SEC ' ) TH(LO=- .3 , SC=. 1 ,NI=6 , ... 

# LI = 1 , UN= ' RAD 1 ) THS ( LO=- . 3 , SC= . 1 , NI=6 , LI=2 , UN= ' RAD 1 ) ... 

# F(LO=- .3 , SC= . 1 ,NI=6 ,LI=4 ,UN= 1 LB 1 ) 

#LABEL (C) TIME HISTORY 

# 

END 

STOP 



72 



X X 



r- 1 m 



LIST OF REFERENCES 



1. Walrath. C.D.. "Adaptive Bearing Friction Compensation Based on Recent 
Knowledge of Dynamic Friction," Automatica, v. 10, n.6, pp. 717-727, 1984. 

Landau, I.D., "A Survey of Model Reference Adaptive Techniques-Thecry and 
Applications" Automatica , v. 10, pp. 353-379, 1974. 

Slotine. J.J.E. "The Robust Control of Robot Manipulators" International 
Journal of Robotic Research , v. 4. n. 2, 1985. 

4. Dynamic Simulation Language/YS, Language Reference Manual, Program 
Number 5798-PXJ (SH20-6288-0), First Edition, IBM Corporation, June 19S4. 

5. Dynamic Simulation Language/VS. Program Description Operation Manual, 
Program Number 5798-PXJ (SH20-6287-0), First Edition, IBM Corporation, 
June 19S4. 



/ j 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22304-6145 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, CA 93943-5002 

3. Naval Weapons Center 1 

Head. Code 394 

China Lake, CA 93555-6001 

4. Naval Weapons Center 1 

Head, Code 3941 

China Lake, CA 93555-6001 

5. Chairman, Mechanical Engineering 4 

Na%'al Postgraduate School 

Code 69 

Monterey, CA 93943-5002 

6. Prof. L. W. Chang, Prof. Code 69CK 1 

Naval Postgraduate School 

Monterey. CA 93943-5002 

7. Mr. Joseph Chan 5 

Naval Weapons Center 

Code 3941 

China Lake, CA 93555-6001 

8. Mr. R. Werneth 1 

N.S.W.C. 

Code U25 

White Oak, MD 20910 



74 



TV'jm Y T~~T‘' s ~ rr T AT"7 



B002 




Thesis 

C3657 Chan 

c.l Effect of friction and 

control parameters on the 
tracking accuracy of a 
target seeker. 




thesC3657 

Effect of friction and control parameter 




3 2768 000 77204 0 
DUDLEY KNOX LIBRARY 




